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Q ■ Abstract 

■ It has been proposed that a sterile neutrino Vh with m^ ~ 50 MeV and a dominant 

decay mode ~^ t^I may be the origin of the experimental anomaly observed at 
^ _ LSND. We define a particular model that could also explain the MiniBooNE excess 

\ consistently with the data at other neutrino experiments (radiative muon capture at 

TRIUMF, T2K, or single photon at NOMAD). The key ingredients are (i) its long 
lifetime {t^ ~ 3—7 x 10"^ s), which introduces a 1/E dependence with the event 
energy, and (ii) its Dirac nature, which implies a photon preferably emitted opposite 
to the beam direction and further reduces the event energy at MiniBooNE. We show 
that these neutrinos are mostly produced through electromagnetic interactions with 
nuclei, and that T2K observations force BR(z^;i — t- VtI) < 0.01 ~ BR(z^/i — t- ^^7). The 
scenario implies then the presence of a second sterile neutrino u^' which is lighter, 
longer lived and less mixed with the standard flavors than Uh. Since such particle 
would be copiously produced in air showers through i^h decays, we comment on 

the possible contamination that its photon-mediated elastic interactions with matter 
could introduce in dark matter experiments. 



1 Introduction 



Although we have today evidence that neutrinos have masses and mixings yp, there are still 
basic questions that past experiments have been unable to answer. The most important 
one concerns their Dirac or Majorana nature: Are neutrino masses purely electroweak (EW) 
or they are revealing a new fundamental scale? Other questions like CP violation or the 
presence of extra neutrino species often use minimality as guiding principle, an argument 
that has worked well in the quark sector. In this sense, the three-flavor framework with 
Am?2 ~ 7.9 xlO"^ eV^, Am^g ^g ^ 2.4x10^3 eV^ and mixings sin^ ^12 ^ 0.30, sin^ ^23 ~ 0.50, 
sin^^ig ^ 0.01 seems to fit well the data from solar, atmospheric and reactor experiments. 

This basic picture, however, has faced a series of persistent anomalies in experiments with 
neutrino beams from particle accelerators. Basically, muon neutrinos of energy below 1 GeV 
seem to experience an excess of charged-current (CC) interactions with an electron in the 
final state. The interpretation of these events in terms of z/^ — )■ oscillations would require 
the presence of (several) sterile neutrinos [2] , and it fits only marginally the combined results 
of LSND, KARMEN and MiniBooNE. 

Here we will discuss a very different possibility. In particular, we will explore a variation 
of the heavy-neutrino model proposed in [3] to explain the LSND/KARMEN anomaly. LSND 
jl] was designed to observe the interaction of muon antineutrinos of E^, < 52 MeV after a 30 
m fiight. The results revealed an excess that was interpreted as z/^ — )■ Ue oscillations followed 
by a CC quasi-elastic interaction giving an observable e"*" and a free neutron. The neutron 
would then be captured to form a deuteron plus an also observable 2.2 MeV photon. The 
LSND signal, however, was not confirmed by KARMEN using a similar technique. 

In [3] Gninenko makes a very compelling case for a massive neutrino as a possible solution 
to the LSND/KARMEN puzzle. In addition to the z/^'s from the decay-at-rest of fi~^ leptons, 
the LSND fiux included 60-200 MeV muon neutrinos from the decay-in-fiight of 7r+ mesons. 
In KARMEN, however, the location of the detector (defining a large angle relative to the 
incident proton beam) eliminates these neutrinos, selecting only antineutrinos of E < 
A 40-80 MeV sterile neutrino would then be above the production threshold there but 
not at LSND. Gninenko shows that if Vh has a sizable component along the muon fiavor 
(|t^/i/iP ~ 10^'^-10~^) it could explain the LSND excess provided that it decays fast enough 
(t/i < 10^^ s) into a light neutrino z/, plus a photon, 

I'h i^a ■ (1) 

He argues that the Cherenkov light from a photon converted into a e^e^ pair would be 
indistinguishable from that of an electron, and shows that the neutron hit by the initial z/^ 
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has the right recoil to provide the correlated 2.2 MeV photon. 

In Gninenko's scenario the production of the heavy neutrino is fixed by the mixing f/^^, 
whereas the lifetime depends on the electromagnetic (EM) coupling /itr (see next Section). 
One can then easily estimate its effects in other experiments. In particular, Uh will be 
produced at MiniBooNE via Z exchange with the nucleons in the detector. Gninenko finds 
[3] that in order to have an impact there its lifetime must be reduced to < 10~^ s, a 
decade below the maximum value suggested by LSND. The large Ufj,h mixing required to 
explain these two experiments, however, will also introduce CC transitions (/i — )■ u^) with 
experimental implications. In |6] it is shown that the model would confiict with observations 
of muon capture with photon emission at TRIUMF [7] , implying a signal well above the 30% 
excess (versus the standard model value) deduced from the data. 

Here we propose the possibility that Uh has a longer lifetime, of order 3-7 x 10~^ s, 
and that its dominant production channel at MiniBooNE is through photon (instead of Z) 
exchange with matter. The main effect of an increased is easy to understand. At TRIUMF 
the target volume is much smaller than cTh, and the number of events with Uh decaying inside 
the detector, 

fi~p ^ ^ j^ijn, (2) 

is proportional to l/rh- Therefore, an increase in the lifetime by a factor of 3-7 respect to 
the value assumed in [6j will reduce the number of events by the same factor. The larger 
lifetime will also reduce the number of events at MiniBooNE. However, as noticed in [6j, the 
EM couplings necessary to mediate decay may also have an impact on its production, 
and we will use them to compensate this reduction. In addition, since the decay length 
of the heavy neutrino grows linear with its energy. 



Ad = CTh- 



(3) 
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the effect of a larger will be stronger on high than on low-energy events at MiniBooNE. The 
energy dependence that we will obtain fits better the data there than neutrino oscillations 
or the prompt zz/^-decay hypothesis by Gninenko. 

In the next sections we define a model for the heavy neutrino and explore its implications 
at MiniBooNE. In section 4 we study the impact of the model on T2K, which further 
constrains the scenario and defines a working model. Finally, in section 5 we sumarize 
our results and discuss the possible implications of the scenario in other experiments. 
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Figure 1: Heavy neutrino production vZ — )■ v^Z through photon exchange with a nucleus 
(left) and decay — )■ (right). 

2 The heavy neutrino 

We will assume that is a Dirac particle of mass rrih = 50 MeV with a left-handed compo- 
nent slightly mixed with the muon neutrino: 

z/^ = cos 6 Vh + sin 6 z/^ , 



where sin^^ = U^h and |f/^/iP = 0.003. We will also assume that the heavy neutrino does not 
introduce any new sources of lepton-number violation. If the light neutrinos z/j {i = e, fi, r) 
are Dirac particles this assumption will not have any implications on the Ui — )■ Uh EM 
transitions. If they are Majorana, however, it implies a relation between the electric and the 
magnetic dipole transitions, 



where we have dropped the prime to indicate mass eigenstates and have taken a real value 
for the coupling fil^. The Lagrangian above is then CP conserving, it only breaks parity. 
We will use it to calculate the production and the decay (in Fig. 1) of the heavy neutrino Uh. 
Since the masses of the light neutrinos are negligible at MiniBooNE energies, lepton number 
will be preserved both in EW processes and in the production or the decay of Uh, and the 
light neutrinos z/j (z7j) will always appear with negative (positive) helicity. 

Production cross section. The production of the heavy neutrino in z/^ collisions with 
matter will include Z- and photon-mediated processes. The first ones, with an amplitude 
proportional to U^h, will couple to the protons and the neutrons in the target oil (CH2) 
at MiniBooNE (collisions with electrons are negligible). We will show, however, that long- 
distance processes mediated by the massless photon, with an amplitude proportional to /i^j., 
may be the dominant ones. 



— sin 6 Uh + cos 9 u. 



(4) 




(5) 
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Let us consider the quasi-elastic production of through photon exchange when a neu- 
trino of energy E^, scatters off a nucleus of charge Z and mass M. The initial neutrino results 
from a meson or a muon decay, and it will always be polarized against the beam direction 
(i.e., it has a negative helicity). If the target is a spin 1/2 nucleus the cross section reads 



daiu^Z 7^h) _ aZ'F\t) (O^ f{s,t) 



(6) 



dt 2 t^{s - M2)2 

where s = + 2EyM, t = {py^ —PuY, F(t) is the form factor, and 

/(s, t) = -2t{s - M^){s -M^ + t)+ mlt{2s + t) - ml{2M^ + t) . (7) 

If the target is a spin nucleus (like ^^C) the cross section can be obtained just by replacing 
the term —m\t in f{s,t) by — t^(t + ml/2). The Mandelstand variable t is simply related to 
the recoil energy T of the nucleus, 

T = (8) 

2M' ^ ' 

whereas the scattering angle 6 in the lab frame is 

fp — T — MX. — 

cose= " . (9) 

+ _ 2EyT - ml 

Notice also that if the neutrino were massless, it would be always produced as a right- 
handed particle of positive helicity, since the EM transition would flip the chirality. The final 
state with negative helicity gives a contribution of order {mh/ E^f to the cross section above. 
For m/j = 50 MeV, at 700 MeV we find that only one out of 5 x 10^ massive neutrinos will 
have the spin against its momentum. In contrast, the produced at MiniBooNE through 
Z exchange will have the opposite (negative) helicity. 

In addition to the helicity of the final neutrino, the main difference between Z and 
photon-mediated processes is the typical distance (or equivalently, the = —t or the recoil 
T) in the colhsion. In Fig. 2-left we plot the distribution {{T/a)da/dT) for u^p — )■ 
at E = 700 MeV. We find that the average recoil is T = 150 MeV when the collision is 
mediated by a Z boson but just 2 MeV in processes that go through photon exchange. 

For form factor F{t) describing the target we have used the expresions in [8]. At 
between 10~^ and 10~^ GeV^, which are dominant in these photon-mediated collisions, the 
incident neutrino interacts coherently with the whole nucleus. At lower values of electrons 
tend to screen the nuclear electric charge, whereas at higher values the collision is mainly with 
the protons inside the nucleus (we use then the usual EM form factors [9]). In Fig. 2-right 
we plot the total cross section for the quasi-elastic process CH2 — )■ I'h CH2 for U^f^ = 0.003 
and fi'^j. = 2.4 x 10"^//^ (see below). 
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Figure 2: Left: for v^p — Vhp through photon (sohd) and Z (dashes) exchange for 

= 0.7 GeV. Right: Total cross section (t(z/^ CH2 ^ Uh CH2) for U^,^ = 0.003 and 

/^i; = 2.4 X io-Vb- 

Decay rate. For the decay — )■ z/j7 (in Fig. 1-right), let us consider a neutrino 
at rest and with spin |+) along the Z direction. It is easy to find [TOl [H] the angular 
distribution of the final photon: 

- ml (1 -cos^) (10) 



dcos6' 327r 

where //^r = Z]i(/^tr)^ 9 is the angle between the photon momentum and the Z axis. 
Therefore, the photon prefers to exit against the spin of the initial heavy neutrino. Anal- 
ogously, when the particle decaying is an antineutrino Uh the final photon will be more 
frequently emitted along the direction of the spin. 

For a neutrino mass rrih = 50 MeV, the lifetime r^, = 5 x 10~^ s is obtained with 
fitr = 7.2 X 10~^ GeV~^ = 2.4 x 10"^^. This is a relatively large value of /itr, as it will be 
generated at one loop and must include a fermion-mass insertion. The ultraviolet completion 
of our Uh scenario would then require the presence of extra physics at the TeV scale (see, for 
example, the left-right symmetric models in |llj). 

When the decay Uh — )■ i^ij includes several neutrino species {i = fi, ...) the branching 
ratios are just BRj = (/ijj,)^/ /i^j,. Notice also that since the beam at MiniBooNE is composed 
of (mostly) muon neutrinos, the production cross section through photon exchange is just 
sensitive to (/i(^r)^- 
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Figure 3: Flux at MiniBooNE in the neutrino mode for 5.58 x 10^'^ POT (left) and in the 
antineutrino mode for 11.27 x 10^'^ POT (right). 

3 Heavy neutrino events at MiniBooNE 

MiniBooNE [12] has run in neutrino [13] and antineutrino [2] modes. We will first analyze 
the results for neutrinos from 6.46 x 10^° protons on target (POT) presented in [13]. These 
data were initially used by MiniBooNE to exclude the neutrino oscillation hypothesis favored 
by LSND (see below), although latter analyses emphasized the anomaly observed at low 
energies [I3] . We will also estimate our prediction for MiniBooNE in the antineutrino mode 
with an exposure to 11.27 x lO^'' POT [15], where the data seems more consistent with the 
3+1 neutrino-oscillation picture. Our fit to the total muon neutrino flux dA^,^/(dS'd£'j,) [12] 
in these two cases is given in Fig. 3. 

MiniBooNE is a sphere of fiducial radius i? 5 m filled up with CH2 of density p = 0.86 
g/cm^. When a neutrino enters the detector it finds an oil column of length L, with L 
taking values between and 2R depending on the point of entrance (see Fig. 4). To define 
an observable Uh event, first z/^ must interact at a distance / < L from the entrance and 
produce the heavy neutrino, and then Uh must decay within a distance L — I, i.e., before 
leaving the fiducial volume. It is easy to see that, while the probability to interact is just 

Pi = , (11) 

mcH2 

the probability to interact and decay within the distance L becomes 

where the decay length A^, defined in Eq. ([3]), depends on the heavy neutrino energy = 
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Figure 4: Geometry of the detector at MiniBooNE. 



Ely — T and its hfetime r^. Since the photon-mediated processes under study introduce 
small recoils, in our estimate we will consider that Uh is produced forward and with energy 
Eh = Ey. Integrating over all the points of entrance into the detector we obtain that the 
energy distribution of heavy neutrinos decaying inside MiniBooNE is 



dEy dSdEy mcH2 Jo \ 2%V3 

where V = inR^/S, mcH2 = (14/6.022) x 10^^^ g and the energy dependence of the total 
cross section a and of is understood. 

The Uh energy distribution dNh/dE^, above must then be translated into a distribution of 
visible energy. As explained in the previous section, most of the neutrinos Uh produced via 
photon exchange have positive helicity, and when they decay — I'i'j the gamma tends to 
exit backwards. If carries E^, Eq. (fTO|) implies that in the lab frame the energy of the final 
photon will be distributed linearly between -Emin = Ey{l — f3)/2 and -Emax = Ey{l + (3)/2, 



with P = Jl — ml/ E^. The distribution of x = E^J E^, is then just 



f{x) = 2-^^^, (14) 



ram) 



with Xmin(max) = -Emin(max) /-E,.. lu Fig. 5 we plot /(x) for Ey = 0.7 GcV and rrih = 50 MeV. 
The energy distribution dN^/dE^ of the photons from heavy- neutrino production and decay 
is then obtained from the integral 

d^-. - f'd.^^(£,/.)M, (15) 



dE^ Jo dEi, X 
where dNh/dEy is evaluated at E,, = E^/x. 

In [13] the results are presented in terms of the energy of the incident neutrino, which is 



reconstructed from the visible energy and the scattering angle of the final electron. In the 
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Figure 5: Fraction of energy taken by 7 in — )■ 1^47 {E,y = 0.7 GeV, nih = 50 MeV). 
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Figure 6: Left: Energy distribution of heavy neutrinos produced in the detector (dashes), 
of heavy neutrinos decaying inside the detector (long dashes), of photons from v^, — > vf^ 
(dots), and of Vh events reconstructed as CC interactions (sohd). Right: Energy distribution 
of Vh events reconstructed as CC interactions (sohd), of events from neutrino osciUations for 
sin^(26') = 0.004 and Am^ = 1 eV^ (long dashes), and excess observed by MiniBooNE in the 
neutrino mode. 



events under study here the final photon will typically exit forward, defining a small angle. 
For example, in average a 400 MeV neutrino would produce through photon exchange a 
with an angle of 7°, and its decay would add 22° to the final photon trajectory. In contrast, 
a ly-mediated process would imply an electron with an average angle of 65° (see [Ml [IT] 
for a study of quasielastic neutrino scattering). In Fig. 6-left we plot the energy distribu- 
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Figure 7: Energy distribution of events reconstructed as CC u interactions (solid), of 
events from neutrino oscillations for sin^(2^) = 0.004 and Am^ = 1 eV^ (long dashes), and 
excess observed by MiniBooNE in the antineutrino mode (with error bars). 

tion dN/dE of neutrinos produced inside the detector, of neutrinos that decay inside 
the detector, of photons from Uh decays, and of the neutrinos that one would reconstruct 
assuming that the visible energy comes from a CC interaction. These results include the 
events where Uh is produced through Z exchange, which for U'^f^ = 0.003 account for just 
a 12% of the total. We have taken Th = 5 x 10~^ s and BR^ = 0.01, which correspond to 
= 2.4 X 10-Vb- 

In Fig. 6-right we summarize our results for MiniBooNE in the neutrino mode. We give 
the energy distribution of the events reconstructed as ly-mediated interactions together 
with the distribution expected from neutrino oscillations for sin^(2^) = 0.004 and Am^ = 
I eV^ and the excess observed at MiniBooNE. These oscillation parameters had been favored 
by the LSND anomaly, but they imply a 500-700 MeV excess that was initially excluded by 
MiniBooNE. In contrast, the long-lived heavy neutrino hypothesis seems consistent with the 
data. 

In Fig. 7 we plot the results from an analogous analysis for 11.27 x 10^*^ POT in the 
antineutrino mode of MiniBooNE [T5]. The data at 500-700 MeV seem to favor the LSND 
oscillation hypothesis, although the excess observed at lower energies would still remain 
unexplained. Our long-lived heavy neutrino, instead, provides a reasonable fit at all energies. 
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4 T2K implications 



LSND and MiniBooNE led us to define a model for Uh where the EM transition that fixes 
its lifetime, = 2.4 x lO^^/i^, is 10 times larger than the coupling required to produce 
it at MiniBooNE, = 2.4 x lO'^/i^. This implies that the branching ratio for the decay 
into the muon-neutrino flavor {uh — ?■ is just around 1%. The obvious question is then 
what other flavor may account for 99% of the decays, and the tau neutrino seems the natural 
candidate. 

To check the consistency of this possibility one needs a I'r beam, and the long-baseline 
experiment T2K [TSl HH] provides such a beam. Muon neutrinos are produced at Tokai with 
El, ^ 0.6 GeV and they travel 295 km to Super-Kamiokande, where most of them have 
oscillated into the tau flavor: 

, ^ . 1.27 Am2 (eV^) L(km) , , 

p.,^., ^ sin^ 2^23 sin^ E%eV) ' ^^^^ 

with sin^ 26'23 ~ 1 aiid Am^s ~ 2.4 x 10"'^ eV^. Therefore, T2K would be able to measure 
the angle through oscillations of muon into electron neutrinos, 

■ 2n .20/, . 2 1-27 AmL(eV^) L(km) 
^ sm2 023 sm2 2^13 sm2 El{Gev] ' ^ ^ 

but also the anomalous EM interactions of Ur with matter to produce the heavy neutrino i^h 
j20] (see also [21] for possible effects of Uh on tau physics). 

The inner detector (ID) at Super-Kamiokande is a cylindrical tank with a fiducial volume 
of radius i? = 17 m and 32 m high {i.e., 22.6 tones of pure water). An outer detector (OD) 
enclosing the tank excludes events with activity at distances below 4.5 m from the fiducial 
volume in the ID. Therefore, there are two ways to generate a signal through decays. The 
heavy neutrino could be produced in the rock before the detector, fly through the OD and 
decay inside the fiducial volume in the ID, or it could both be produced and decay inside 
the ID. 

If only muon neutrinos can produce Uh through photon exchange [i. e., fi^,. = 2.4 x 10~^fiB, 
= 0), the energy distribution of the heavy neutrinos that are produced and decay inside 
the ID is 




dEi, dS dEj, fTiH^o Jo ^ \ 2Rcos9 

where is defined in Eq. (jS]), o" is the cross section to produce a Uh in the collision of z/^ 
with water, p = 1 g/cm^, and we have taken the neutrino flux dN^y/dS dEi, in |19j. For 
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comparison, the distribution of v^, producing electrons in collisions with water would just be 

(19) 



The calculation of the energy distribution of the neutrinos produced in the surrounding 
rock (beyond a distance Dq = 4.5 m of the ID) and decaying inside the detector is also 
straightforward, 

dE^ dSdE^ rrir n R Jo \ J ^ ' 

where we have taken = 3.2 g/cm^. 

We obtain that whereas z/^ — )■ oscillations with sin^ 26*13 = 0.1 would imply 6 events 
of energy above 100 MeV, i/^ colhsions producing a Uh that decays into z/j7 would introduce 
1.1 events, 75% of them from heavy neutrinos created in the rock and decaying inside the 
detector. These long-flying neutrinos dominate over the ones produced inside the detector 
due to the larger density of the soil sorrounding the detector and the relatively long lifetime 
{Xd = 15 m for Eh = 500 MeV) of Uh- Therefore, in our framework the MiniBooNE anomaly 
could also contribute to the signal observed at T2K [19], implying a smaller value of the 
mixing ^13. The contribution would be larger for a non-zero value of of the coupling fil^. 
that defines the Ur — ?■ i^h EM transitions. Since these Uh events have different kinematical 
and energy distributions, it seems clear that an increased statistics there could disregard 
this possibility. Incidentally, the events with Uh produced in the rock would be distributed 
preferably near the point of entrance into the detector (the distribution is proportional to 
g-o/Ad^^ as it has been observed, whereas the — )■ z/g oscillation hypothesis implies a flat 
distribution. 

The previous result also implies that /i^j, < /ifj.- In particular, if the channel Uh — ?■ 
accounted for 99% of the Uh decays, then the number of electron-like events at T2K from 
ly^-Z — 7- UhZ EM transitions would be around 240, well above the 6 events that were observed. 
As a consequence, to be viable the model requires an additional sterilj^ neutrino which 
is lighter than Uh and has the coupling /x^j.' ^ 2.4 x lO'^/i^ required to mediate Uh — )■ Uh'^y at 
the right rate (r/j ^ 5 x 10~^ s). 

Finally, we would like to comment on the signal that the model could imply at the near 
detector (ND280) [18] in T2K. When the neutrino beam reaches the off-axis detector at 
ND280 it has not changed its muon flavor yet. A typical event would consist of an initial 



^The electron flavor seems also disfavored, since the process VeZ — > lyhZ with Vh — > i^ej in the atmosphere 
would introduce a signal identical to CC interactions. 
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interaction producing the heavy neutrino and a small energy deposition (a 10 keV-10 
MeV nuclear recoil) in the Pi-zero detector, followed by the single photon from its decay at 
any point in the tracking system (three TPCs and two thinner FGDs). We estimate around 
2.1 events of this type per 1000 CC interactions, plus 0.8 events with the Vh both being 
produced and decaying inside the Pi-zero detector. In the tracker the photon would convert 
into an e^e^ pair that could be distinguished from the single electron plus recoil in a CC 
interaction, whereas events with neutral pions giving two photons seem also clearly different. 
Therefore, ND280 may provide some ground to probe the MiniBooNE/LSND anomaly and 
establish its electron or photon origin (if any). 

5 Summary and discussion 

LSND observed an excess of 3 x 10~^ interactions with an electron in the final state per 
each CC event. The beam in this experiment included antineutrinos of £^ < 60 MeV 
from the decay-at-rest of ji^ leptons but also neutrinos of up to 200 MeV. KARMEN tried 
then to confirm the LSND anomaly using a beam that excluded the high-energy region of 
the spectrum, and it did not see such an excess. 

Gninenko has proposed that the LSND anomaly could be explained by a 50 MeV neutrino, 
a mass which is beyond the reach at KARMEN. This hypothesis requires that Vh decays into 
Vi^ with T/j < 10~^ s and that the production mechanism is mediated by the Z boson (or 
other massive particle), as the in the collision must be large enough to produce the free 
neutron also detected at LSND. 

A similar rate of anomalous interactions (3 x 10~^ per t/^ CC event) has also been observed 
at MiniBooNE both in the neutrino and the antineutrino modes. In contrast with LSND, the 
energy distribution of these events does not seem consistent with a simple 3-1-1 oscillation 
scheme, as it peaks at low energy (where the i/^ — > probability vanishes) and is non- 
significant at the expected oscillation maximum. The initial neutrino data were actually 
used by MiniBooNE to exclude the oscillation hypothesis favored by LSND. Although the 
first MiniBooNE results in the antineutrino mode were different and could favor a 3-1-1 
scheme, the increased statistics obtained during the 2011 run (in Fig. 7) shows consistency 
with the observations in the neutrino mode (in Fig. 6-right). 

Gninenko 's hypothesis can have an impact at MiniBooNE only if the lifetime of the 
heavy neutrino is reduced below 10~^ s (versus the maximum 10~^ s required at LSND). 
The energy distribution that it implies is better than the one obtained from oscillations, but 
probably still a bit flatter than the one observed at MiniBooNE. In Fig. 8 we plot Gninenko's 
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Figure 8: Energy distribution of Vh events in the Gninenko model (|f/^hp = 0.007, = 10^^ 
s, fi^^ = 0) (dashes), in our set up {\U^h\^ = 0.003, Th = 5x 10~^ s, fi^v = 2-4 x lO'^/is) (sohd), 
distribution of events from neutrino oscillations for sin (29) = 0.004 and Am^ = 1 eV^ (long 
dashes), and excess observed by MiniBooNE in the neutrino mode. 

model for a Majorana Uh with |f//i/ip = 0.007 and Th = 10^^ s. If Uh were a Dirac particle its 
distribution would be even flatter, as i>h has then negative helicity and decays predominantly 
into a forward photon. At any rate, the mixing and the lifetime in Gninenko's model implies 
too many events with radiative muon capture at TRIUMF [6]. 

Our initial observation is that the coupling required to explain the decay of may also 
imply its production through photon exchange with matter. While the lifetime is a measure 
of /i^j, = J^ii^J'trYy where i runs over all the light neutrino flavors, the production rate at 
MiniBooNE depends on fi^j. only. The photon-mediated production of Uh will not contribute 
significantly to the anomaly at LSND, as the interaction is unable to free a neutron from its 
nucleus (it involves charged particles and is very soft). 

Therefore, we keep the mixing If/^^hp ~ 0.003, required to explain LSND through Z- 
mediated interactions, but we also increase the lifetime by a factor of 3-7. First of all, 
this reduces in the same proportion the number of CC events with decaying inside the 
detector at TRIUMF, providing for consistency with this experiment. Second, the number 
of Uh decays at MiniBooNE will also be reduced, but this can be compensated by the larger 
number of heavy neutrinos obtained from photon-mediated processes. Finally, the longer 
lifetime implies that low-energy neutrinos are more likely to decay inside the MiniBooNE 
detector than the more energetic ones. In addition, the negative helicity of the initial z/^ 
implies that the EM transition produces a (Dirac) polarized forward, with positive helicity. 
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and that the final photon is emitted preferably backwards. These two factors define a 
spectrum that seems consistent with the observations at MiniBooNE. 

The scenario implies that Uh — ?■ accounts for just around 1% of the decays, and T2K 
forces that the branching ratio for Uh — ?■ is even smaller. In any case, an increased 
statistics at T2K could be used to probe the model. Other experiments could also put 
important constraints. Most notably, NOMAD [22] has recently set bounds (not discussed 
in previous sections) on single photon events in neutrino interactions at ^ 25 GeV. Although 
Gninenko [23] has shown the consistency of his model with the data, our scenario involves 
larger lifetimes and production cross sections for the heavy heavy neutrino. The increase in 
Th will introduce two competing effects that tend to cancell each other: the heavy neutrinos 
produced in the rock are less likely to decay inside the detector, but they can be produced 
further from the detector and still reach it. We estimate that for Th = 5 x 10~^ s our model 
implies around 2 x 10~^ events per each i/^ CC event (above the 1.6 x 10~^ limit established 
there). While lower values ^ 3 x 10^^ s would be safer, a more definite statement would 
require a full MonteCarlo simulation. 

Since the dominant decay channel — )■ i/j7 is not into the muon nor the tau fiavors 
and the electron fiavor seems also disfavoured, an interesting possibility would be that the 
heavy neutrino decays into a photon plus another sterile heavy neutrino Uh'. The mixing of 
this second neutrino with the standard fiavors could be much smaller than the one required 
for Uh (|f//x/iP = 0.003), while a 5-10 MeV mass and a lifetime below 1 s would make it 
consistent with bounds from astrophysics and cosmology [2l]. It could decay Uh' — ^ 'l^e,fi,T 
with cTfii > 10 km. Such a long lifetime would make it very frequent in the atmosphere 
[25] : it would be produced in Uh ^h'l, with the parent Uh appearing in ^ 0.3% of all 
kaon and muon decays [3]. Such neutrino could also have a relatively large EM dipole 
moment, able to mediate elastic interactions u^'Z — )■ Uh'Z with a cross section softer but 
orders of magnitude larger than the ones mediated by the Z boson. In particular, at 1-100 
MeV energies one could expect collisions with the typical recoils (few keV) in dark matter 
experiments. We think that it would be interesting to analyze to what extent such neutrino 
could contaminate experiments like DAMA/LIBRA [26] or CoGeNT [27]. A similar effect 
caused by solar neutrinos oscillating into a sterile flavor has been discussed in [28l [29] 

To conclude, although the current framework for neutrino masses and mixings explains 
most of the experimental data, the picture that we have today is still far from complete. 
The possibility that neutrinos involve new physics at the MeV, a scale that we have been 
exploring for many decades, may sound unlikely or unappealing. We think, however, that it 
should be considered, specially as long as the current experimental anomalies persist. 
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